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I. INTRODUCTION 



Quantum field theories (QFT) of interest to particle physicists are plagued by UV di- 
vergences. Typically, we are confronted with theories which favour local interactions, which 
means that we can remove these divergences order by order in the loop expansion. This iter- 
ative removal of divergences, all the while maintaining the locality of the theory, has to fulfill 
certain combinatorial properties, succinctly summarized by Zimmermann's forest formula. 
It was a great achievement when the self-consistency of the whole renormalization procedure 
was proven While there is little doubt that renormalization theory allows one to obtain 
sensible answers from an a priori ill-defined theory, only recently has it become clear that 
its underlying mathematical structure is by no means accidental, for it relates QFT to basic 
mathematical notions, well-known in low dimensional topology, number theory and 
even in numerical analysis 0. 

In this paper, we want to show to what extent the renormalization of Feynman graphs 
can be incorporated in the language of generalized iterated integrals, and how the intricacies 
of full-fledged QFT, including spin and Lorentz structure, still can be described in terms of 
such generalized iterated integrals in the spirit of pi. 

We will begin by describing how the Hopf algebra structure, associated with renormal- 
ization, can be related to rooted trees, through the nature of the Feynman graphs. In that 
opening section we also show how two kinds of renormalization schemes can be connected by 
a convolution product; this convolution is the group product naturally related to the Hopf 
algebra, and dutifully delivers the renormalization group. The next section explains how 
rooted trees acquire decorations, according to the topology of the Feynman diagram, and 
how their evaluation is completely determined by vertex weights and generalized tree fac- 
torials. The following section concerns the representation of Feynman diagrams as iterated 
one-dimensional integrals and how this can be tied to the renormalization group, through 
appropriate choices of end-points. Finally we treat the case of massless Yukawa theory and 
quantum electrodynamics in four dimensions; the former sets the scene for tackling scalar 
theories, while the latter contains the necessary matrix generalization to spin, where the oc- 
currence of form-factors can complicate the argument. The paper ends with a brief section 
embodying our conclusions and highlighting topics for further development. 

II. FORMULARY FOR RENORMALIZATION 

In this section we will briefly describe the Hopf algebra structure of renormalization. We 
start with a basic formulation directly on graphs and then introduce the correspondence to 
rooted trees. 

A. Basics 

Let us consider a Feynman graph T as a set, consisting of vertices and edges of several 
sorts. The edges correspond to various types of propagators, representing inverse differen- 
tial operators of free relativistic wave equations with boundary conditions in accord with 
causality. The vertices correspond to local interactions associated with a Lagrangian which 
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is a polynomial in fields and derivatives. In this way, at the same time T denotes a graph 
which we can draw on paper as well as a unique analytical expression. 

Regarding T as a set of edges and vertices, we can formally consider the power set of 
this set, and distinguish those elements of this power set which themselves correspond to 
superficially divergent Feynman graphs. 

Such a Feynman graph, together with all its superficially divergent subgraphs, constitutes 
a Hopf algebra structure which is isomorphic to a Hopf algebra of decorated rooted trees. 
This isomorphism is described in Fig. [I] and is detailed in 0,|J. (For the reader's convenience, 
we have added an Appendix which summarizes the crucial points.) By this isomorphism, 
the standard Hopf algebra operations on rooted trees correspond to established notions of 
physicists. 

Admissible cuts on the sum of rooted trees representing T are in one-to-one correspon- 
dence with divergent subgraphs. Those are determined by power-counting, thereby allowing 
detection of all sectors in the analytic expression T which correspond to divergent subinte- 
grations in their own right, while the forests of renormalization theory then correspond to 
arbitrary cuts [@-§J]. 

This viewpoint leads to a Hopf algebra structure 0-0 , which we can succinctly formulate 
on overall divergent Feynman graphs T as 



A[1] = 1®1 (1) 

A[r]=£7®r/ 7 (2) 

S[l] = 1 (3) 

S[T] = -r - £ 7 ® r/ 7 (4) 

= -r- £(-i) n ^F®r/F (5) 

e[l] = 1 (6) 

m = o. (7) 



In this notation, the sum XL c r is a sum over all superficially divergent subgraphs, including 
the case 7 = T and 7 = 0. r/7 is obtained by the contraction of 7 C T to a point in T. 
The algebra element corresponding to the empty set is, as usual, the unit 1 of the algebra. 
Hence, the coproduct determines all sectors in the graph which are ill-defined and require 
renormalization. The sum X 7 c r * s a similar sum, with the two cases 7 = T and 7 = 
excluded. 

T is the set of all forests of a Feynman graph, and this set is in one-to-one correspondence 
with all cuts at the rooted tree t-p representing the graph T, while divergent subgraphs are in 
one-to-one correspondence with all admissible cuts, as already mentioned, rip is the number 
of boxes of the forest F, or the number of elementary cuts under the above correspondence [|J . 
Since formulae (l)-(8) are rather abstract for the uninitiated reader, we give an illuminating 
example in Fig. [l|; a proper derivation of these correspondences can be found in 0-§f while 
the Appendix summarizes the crucial notions. 

Hopf algebras, as embodied in (l)-(8) are Hopf algebras of rooted trees for an appropri- 
ately chosen set of decorations With this Hopf algebra structure, the naive renormalized 
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function T[T] associated with T vanishes: 

T[r] := m[(S <g> id)A[T\] = e[T] = 0. 



(8) 



Spelling the formula out, one indeed finds that each term appears twice with opposite signs. 
This must be so as the above formula employs the graph itself as the counterterm of a 
divergent graph, and hence eliminates all divergent graphs. Equation (8) represents a well- 
defined renormalization scheme, but is not of much use in practice. It has a few merits in 
the large Nc expansion however M. 

Usually we want to subtract in accordance with certain renormalization conditions, and 
hence our counterterms should surely remove divergences of integrals but may not completely 
nullify all superficially divergent bare Green functions. We have to choose non-trivial renor- 
malization schemes to obtain non-vanishing but finite results. The essential modification is 
to allow for a map R which modifies the bare unrenormalized Feynman graph without chang- 
ing its divergence. In |J it was shown that all renormalization schemes can be described on 
the same footing, allowing for tree-dependent variations of dimensionful parameters in bare 
Green functions. 



B. Schemes and scale transformations 



We can use the results of || to recast the behaviour under a change of parameters and 
the dependence on renormalization schemes in a nice manner. We first modify the antipode 
via a renormalization map R, and in consequence get a finite renormalized Green function 
T R [T] The R dependent antipode, the Z-factor S R [T] of a bare expression T and the 

renormalized Green function r^[r] explicitly read 



S R [T) = R 



R 



r - E s*[7]r/7 



E (-l) nF+1 F R ( lF )T/ lF 
T R [T] =m[(S R ®id)A[T]\ 

= [id-R] (j2(-l) nF F R (l F )T 



hi 



(9) 

(10) 

(11) 
(12) 



Here, 7^ is the subgraph corresponding to the forest F, and F R is the evaluation of this 
subgraph using R to evaluate boxes of forests corresponding to cuts, cf. Fig. |[ 
Introducing the convolution product (</> * ip)(T) = m[(0 ® tp)A(T)] we find 



r R [r) = (s R *id)[r], 



(13) 



which makes the renormalized Green function transparent as a ratio of a bare Green func- 
tion T = id(T) to a counterterm S R (T). Hence our notation emphasizes the fact that the 
renormalized Green function is a function of the bare Green function as well as the chosen 
renormalization scheme. A detailed discussion of all its properties can be found in H. Here 
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we only remark that the group structure embodied in the renormalization group is a direct 
consequence of the structure of a renormalized Green function as a convolution product []5],|| . 

Of course, bare functions T are parametrized by masses, external momenta and regu- 
larization (cut-off, dimensional-regularization scale /x, etc.). We can vary such parameters. 
Essentially, the choice of a renormalization scheme amounts to a choice of conditions for ex- 
ternal parameters. Hence, a renormalized Feynman graph is determined by the choice of two 
sets of such parameters, one for the bare diagrams, one for the counterterm diagrams. If we 
want to highlight the dependence of the renormalized Green function on such parameter sets 
i, k say, we explicitly write r^*.; this assigns to a Feynman graph a renormalized expression 
with a specified set of dimensionful parameters k used in the bare Green functions and a 
set of parameters % used for counterterms. For example, an on-shell renormalized off-shell 
Green function will be calculated by using bare Feynman diagrams with off-shell momenta 
and imposing on-shell conditions to fix them in the counterterms. 

Now, from |J, we know that the most general renormalization scheme, including for 
example minimal subtraction (MS), can be captured by suitable modifications of these scales. 
There one discovers that the change of schemes and scales can be succinctly summarized by 
the following formula 

r Ri , k (r) = [T Rid *T R . )k }(r). (u) 

In r^.j, the first subscript at the renormalization map R indicates the choice of scales 
for the counterterms, while the second one indicates the set of external scales used in the 
bare diagrams. Equation (TTJ]) shows how an intermediate set of scales j allows for group- 
like transformation laws. From this, one obtains renormalization group equations, operator 
product expansions and cohomological properties of renormalization || with ease. Depend- 
ing on the choice of parameters to be varied, one arrives at Wilson's viewpoint by varying a 
cut-off || or at the Callan-Symanzik equation from a variation of the renormalization point. 
Details will be given in future work. 



C. Example 



The following example, of a three-loop fermion self-energy contribution in QED in four 
dimensions, will make the above ideas more transparent and their implementation more 
concrete. We first calculate the coproduct for the diagram: 



A 



1 + 1 

+ "4(g) -X^X, ; 

The antipode Sr[ ^<^ k ] is then found to be 



s R [-^) 



-R 
+R 



R 



+ R 


R 


< 




+ R 


R 


< 





R 



R 



R 



< 



Whenever we see R, this means we have to evaluate the analytic expression corresponding to 
the Feynman graph on which R acts according to some chosen rules. For example, in minimal 
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subtraction we would introduce a regularization, say the dimensional method, evaluate the 
bare Feynman graph as some Laurent series in the corresponding regularization parameter, 
and then read R as the projection onto the pole term. On the other hand, in on-shell BPHZ 
renormalization, we would introduce no regularization, but let R modify the integrand in a 
way such that it fulfills certain constraints when external momenta are on-shell. From the 
results in we know that all such schemes are fully equivalent. In particular the result of 
any scheme can be recast into BPHZ form with appropriately chosen constraints for external 
parameters. 

Whatever scheme is chosen, the renormalized Green function for the example finally 
becomes 



Ri 




Here, id is the identity map which leaves the bare diagram unchanged. The same result is 
obtained when one works with rooted trees. 



III. FEYNMAN DIAGRAMS AS REPRESENTATIONS OF ROOTED TREES 

Particularly simple classes of Feynman graphs are given by graphs which are represented 
by rooted trees which have the same decoration at each vertex, and which give a multiplica- 
tive representation defined below. We will soon see that the general case can be reduced 
to similar representations, possibly allowing for appropriate matrix structures to incorpo- 
rate necessary modifications due to spin. Hence it is legitimate to start the study of the 
renormalization problem on graphs as simple as the ones shown in Fig. [| These simplest 
graphs are completely governed by the combinatorics of undecorated rooted trees. Indeed, 
tree-factorials and vertex weights to be defined below are sufficient to determine such graphs 
completely @,||. 



A. Undecorated trees 

Let us first recapitulate some basic notions of rooted trees, following ||. For further 
details we refer the reader to the Appendix. Let T' ' be the set of vertices of the rooted tree 
T with n vertices. For any vertex v of T which is not the root, let t v = P C (T), where c is 
the single cut which only removes the one edge incoming to v (all edges are oriented away 
from the root). If v is the root r, set t r = T. Also, w(v) = #(t v ) is the number of vertices 
of the tree t v . Then, we define the tree factorial by 

T ! := J] (15) 
Hence, for the root, w(r) = n = #(T). Fig. |3| gives instructive examples. 
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In the following, a variable z plays the role of a q 2 /fJ? in field theory, and x plays the 
role of the regularization parameter, and can be identified with (4 — D)/2 in dimensional 
regular ization. 

We now consider representations of Tt defined as follows. Assume that we are given a set 
of functions Bk(x) which are Laurent series in the regularization parameter x, each having 
a first-order pole. Using vertex weights as defined above we can construct a function 

Mt)-= II B w(v) {x)z~ nx = B t {x)z~ nx , (16) 

where n is the number of vertices of t and z is to be regarded as the scale parametrizing 
the representation. We have also abbreviated rLeTM Bw(v)( x ) — B t (x). Quite a number of 
interesting applications can be brought into this form |5||J. Typically, one finds such repre- 
sentations whenever one iterates a (scalar) Feynman diagram in terms of itself, as described 
by a rooted tree @,|]. Even if one only iterates one- loop integrals, such representations offer 
a variety of interesting number-theoretic properties realized in the same manner by different 
quantum field theories, ranging from the absence of transcendentals in counterterms derived 
from ladder graphs ]l3|-0 (which represent trees of the form B+(e), ||) to the absence of 



non-rational coefficients in quenched QED [IB] to similar phenomena when summing up all 
diagrams with Connes Moscovici weights ||. 

Thus bare Green functions are given by <p z (t). Then we define the MS renormalization 
scheme by setting 

Rms[<Pz] =< <Pi >, (17) 

where angle brackets denote projection onto the pole part of the Laurent series in x inside 
the brackets, and we set z — 1 inside these brackets. Next, let the minimal subtracted 
counterterm be defined by 

s Rli3 (<j>,)(t) = -< <p z (t) + E s RMS (Mki))Mk2)) >, (is) 

which is the direct translation of (|9p using the map from Feynman graphs to decorated rooted 
trees, and let the renormalized Green function be defined by r M>s (0 2 )(t) = [SR MS o(p z ic<p z ](t), 
in accordance with (|i~T|), as expected. The reader should have no difficulties confirming that 
for t 2 , the rooted tree with two vertices, 



S RMS (<j> z )(t 2 ) = -< B 2 B 1 > + «B l >B 1 >:= Z^ s , (19) 



?MS 

and 

rjw.(*a) = B 2 B lZ - 2x - <B 1 > B lZ - x - < B 2 B X > + « B x > B x > . (20) 

Therefore the antipode Zf' IS := Sji MS {(j) z ){t) in MS reads, for the first few trees, 

Z?1 S = -<B 1 >, (21) 
z ms = _ K B2Bi > + <<Bl >B 1 >, (22) 

K[ S = -< B^Bi > + «Bi> B 2 B X > + « B 2 B X > B x > 



«< B x > B x > B x >, (23) 
< B Z B\ > +2 « B x > B 2 B X > 

«< B x > B x > B x >, etc. (24) 
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which can be summarized by the formula |J 



E 



< 



full cuts C of t 



n<^ 



> 



(25) 



This being the result in the MS scheme, it was shown in || how to express it in terms of 
a BPHZ type scheme, where the use of a regularization scheme can be completely avoided. 
Translated to the above representations, what we want to achieve is the following: the above 
functions representing rooted trees depend on a function B t (x), determined by a chosen tree, 
a regularization parameter and the chosen one-loop graph (the decoration) which is iterated 
according to the tree structure. Further, there is the scale z, delivered in simple bare 
Feynman diagrams typically by a variable, such as the square of a single momentum, q 2 . 



An on-shell scheme would simply subtract at a chosen point q 2 / fi 2 



1, for example. 



An MS scheme would also modify the Laurent series B t , by projecting to < B t >. In the 
former case, one typically obtains the difference of integrands evaluated at scales z versus 
1, and hence one may subtract integrands before having to regularize. In the latter case, 
one modifies Laurent series after regularization and integration. We want to show how 
appropriate subtraction points fit can be defined so that the MS renormalized functions equal 
on-shell renormalized functions subtracted at fi t - 

The idea is to use tree indexed parameters fi t . The antipode Zf in a subtraction scheme 
using tree-indexed parameter sets fi t reads 



E 



-i 



full cuts C of t 



1 % 



c 



#{t R )x 



(26) 



where P c {t) = Ili which one immediately obtains by standard formulas 0. The first two 
trees ti,tz deliver immediately 



-B t2 fXt 2 2x + B tl B tl fit 2x , 



(27) 
(28) 



as the reader should confirm using the equation (pq) . For the definition of full cuts, appearing 
in the summation above, we refer the reader to the Appendix. 



Equating Z t MS to Z% determines \i t recursively 



fit = exp 



where 



B' :-- 



E (-i) 



no 



full cuts C of t 



C\ 



S RMS {<Pz){tl 



(29) 



(30) 



One also confirms that the MS-renormalized Green function can now be expressed 
through the on-shell subtracted one, where the on-shell condition for <p z (t) signifies sub- 
traction at z = fi t in this case. Thus 

rW0*)(i) = r Ril (,/g (t) (31) 

holds. Also, one immediately derives that B' t /B t = 1 + 0(x), so that we properly arrive 
at well-defined conditions on the external parameters fit to express the renormalized Green 
function in MS through a BPHZ subtracted Green function with tree-indexed boundaries. 
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B. Decorated trees 



What is the situation that confronts us in general? In the previous section, we dressed 
every vertex of a rooted tree with a decoration B w ( v \, and then assigned the function 4> z {t) = 
z -#(t)x y\ B w r v \ to a rooted tree, the product being taken over all vertices of the tree. (Note 
that = n = w(r), where w(r) is the vertex-weight of the root.) Clearly, the vertex 
weight w(r) of a root of an undecorated tree is simply the number n of vertices of the tree. 
This equality does not survive in the general case. 

But if we only very slightly modify this equation and write 

<t> z {t) = z- w{r)x I[B w{v) . (32) 

V 

we obtain the form which permits all necessary generalizations. Such generalizations include 
the following cases: 

• In general, decorations can be provided by analytic expressions of any loop order. We 
will have to attach appropriate weights to the vertices themselves, and hence will have 
to generalize the concept of vertex weights. 

The crucial step is to devise some sensible notation. An arbitrary rooted tree has 
vertices decorated by graphs corresponding to analytic expressions without subdiver- 
gences which themselves can have any loop order. We still assume that these analytic 
expressions can be described by scalar functions. Loop integrations of decorations 
at vertices below the vertex considered will alter the measure by an amount (k 2 ) w ^ 
where V{ is the generalized vertex weight defined in Fig. |4] and k is an appropriate 
external momentum for the decoration at At each vertex we indicate the loop 
number of its decoration. The generalized vertex weight is then the sum of all these 
loop numbers, assigned to the vertex in question and all vertices below it. /,From this, 
one obtains the generalized tree factorial, which is to be used in (|32|). 

• We might consider graphs which do not come from scalar theories. In order to cope 
with form factors, one needs to devise an appropriate matrix calculus; this will be 
worked out in a detailed example in a section below. Essentially, the master equation 
( p2|) remains valid, but becomes now a matrix equation. This is the most non-trivial 
generalization, and we will discuss it in detail presently. 

In the same manner, one can use matrices to keep track of internal structure of a 
graph: at some stage one has to bookkeep the location of sub divergences in a larger 
graph, and this leads to a matrix calculus which is a straightforward generalization 
of the example which we consider below. Indeed, one can enumerate all vertices of 
a graph and can construct suitable matrices which store the information about the 
nature of these vertices and the propagators connecting them. In this manner, the 
previously discussed changes of measures due to sub divergences operate on entries of 
such matrices, and it is a straightforward notational exercise to set-up a convenient 
calculus. We will not consider the most general case here, but rather present below an 
example for QED which should make the idea sufficiently clear. 
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• In general, divergent subgraphs depend on masses and more than one external mo- 
mentum; hence on multiple scales. For example, a vertex subdivergence will depend 
on more than one internal momentum usually and this can spoil the multiplicative 
form of fl32l). This difficulty can be solved by a proper decomposition into sums of 
expressions, each of the form fl32|). This was described already in [[J. One proceeds 
by decomposing any function of multiple scales into one of a single scale and another 
finite multiple scale function. This is always possible by adding zero in appropriate 
ways. 

The final result is that from a local QFT we essentially get analytic expressions which 
can be expressed in the form fl32|). We now utilize this fact to show that they essentially 
behave as generalized iterated integrals. 



IV. FEYNMAN DIAGRAMS FROM GENERALIZED ITERATED INTEGRALS 

Now that we have established that (j32|) is an acceptable prototype expression into which 
Feynman diagrams can be resolved, we want to utilize the fact that the simplest represen- 
tation of such form can be obtained from generalized iterated integrals. Such generalized 
iterated integrals were prominent in || and give a convenient means to investigate renor- 
malization. Generalized iterated integrals are distinguished from ordinary iterated integrals 
by the fact that boundaries of the domains of integration are determined as a function of the 
rooted tree which the integral represents. In this respect, ordinary iterated integrals use a 
function which is constant on rooted trees. For such a generalization, the shuffle symmetry 
of iterated integrals is lost. 

The extra freedom gained by using this generalization will allow us to rewrite Feynman 
diagrams as such generalized iterated integrals. We will see that the generalization only 
affects non-leading UV divergences, and hence we will obtain later the result that the leading 
UV divergences of Feynman integrals still obey a shuffle identity. 

At this point, it might be worth stressing that if we were able to reduce Feynman 
diagrams completely to ordinary iterated integrals we would obtain considerable progress in 
understanding the analytical and transcendental structure of Green functions, as the theory 
of ordinary iterated integrals is quite rich and highly developped. We hope that a similar 
theory will be emerging in the future for generalized iterated integrals. 



A. Iterated Integrals 

Essentially, we want to go one step further as in || and express contributions to Green 
functions in QFT as iterated integrals. We will start our considerations by reminding our- 



selves of some basic properties of iterated integrals 110,11 , specializing to the case of a 



single function f(x) with associated one-form f(y)dy on the real line. We assume that f(y) 
behaves ~ 1/y for large y. 

Iterated integrals built with the help of / are parametrized by an integer n and two real 
numbers a, b say. They are defined by 
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a, ft 



, Va, b G R, 
Hence we can write them as an integral over the simplex, 



7 a!b = / f( X l) ■ ■ ■ f{Xn)dX! . . . dx n . 

J a<xi<...<x n <b 



(33) 
(34) 

(35) 



We will also consider 



a, b 



where T is a rooted tree with B_{T) = ITjtj. 

A generalization which turns out to be quite useful in practice is to let even the boundaries 
a, b be indexed by decorated rooted trees. This is similar to what we did before when we 
expressed the results of an MS calculation by an on-shell scheme with suitably adjusted 
tree-dependent parameters. Here, parameters are provided by boundaries in the iterated 
integral. Hence, we define an iterated integral which serves as a bare Green function 



Gb t ,oo(t) 



bt 



f(x)Y[G Xt00 (ti)dx. 



(36) 



The product is again over all branches U of B-(t), and we set Gft ej0 o(e) = 1. This is the same 
as before, except that we now label the lower outer boundary by a decorated rooted tree t 
and set the upper boundary to infinity. We assume that a renormalization map R a shifts 
the outer lower boundary b t to another value a t , and that the coproduct action extends to 
this label. We define 0& to be the map which sends t — > Gb t)00 (t). 

It is worthwile citing here an example taken from which exhibits the intricate nature 
of a change of scales in full generality. For the rooted tree with two vertices, £2, having 
coproduct 



A(* 2 ) = t 2 ®e + e®t 2 + t l ®t l 
we formally obtain a counterterm 

SRa,{<l>b){h) = -RaiMh) + m[(S Ra o cj) b a 
(where A'(T) = A(T) - [1 <g> T] - [T ® 1]) as 



(37) 



)A'(t 2 )]] 



f{x)f{y)dydx. 



(38) 



and a renormalized iterated integral r^^) = m [{SR a ° <pb <E> 0&)A(t2)] as 



r a ,&(*2) 



00 roo poo poo 

/ + 



f(x)f(y)dydx. 



(39) 



In this notation, a, b are to be regarded as representing actually a whole set of constants 
at,bt, parametrizing the relevant scales for the decorated tree under consideration. 
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We now want to change the set of scales and reexpress how T a ^ can be obtained from 
renormalized Green functions r a s and F S;b . Obviously, the convolution 



r„, fe (t) = [r a , a *r a , 6 ](f) 



holds. This form of Chen's Lemma describes what happens if we change the renormalization 
point. It works for the simple case of iterated integrals in the same manner as for full Green 
functions of QFT. T a ^ b = T a ^ s -k T s ^ b now becomes 



oo roo 



bt 9 Jx 



OO /"OO 



oo roo 



oo roo 



at, Jat 



oo roo 



St. Ja t 



+ 



at 2 Jx 
oo /-oo 



at Jx 



+ 
+ 



oo roo 



oo roo 



oo roo 



oo roo 



f(x)f(y)dydx 
f(x)f(y)dydx 
f{x)f(y)dydx 



roo roo 



f{y)dy 



6*1 



s *i 



f(x)dx. 



This is evidently true, as the reader may readily check. 

Finiteness of T a b now imposes conditions on the tree-indexed parameters, a fact which 
we will utilize in the next section. Actually, we will determine these parameters in a way 
which reproduces the results of a renormalized QFT. 



V. COUNTERTERMS AND RENORMALIZED GREEN FUNCTIONS AS 

ITERATED INTEGRALS 

We are now in the fortunate position of being able to transform Feynman diagrams to 
generalized iterated integrals. To that end, let us consider the simplest case of an iterated 
integral based on the single one-form f(y)dy = y~ 1 ~ x dy, which diverges logarithmically 
when integrated to infinity and x — > 0. 

Then, we obtain the bare iterated integral, 



Gb,oo(t) 



1 



-#(*)* 



t ] -x-#® * 

and associated counterterms Sr ((fib) (t) , which in this case are simply 

S R ^(t) = -R,Mt) + m[(S R ^ ® id)(<fi b ® <f> b ))A'}\ = MS(t)), 
with R^((fib) = <fin- Thus we arrive at the renormalized integral 

r,, b (t) = <fi,( pC (t))MR c (t))- 

Explicitly, 

1 



s Rfl (t 2 ) 



2x 2 



1 ^ 2 

x 1 
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and 

which is nonsingular when x — *■ 0. 

Next determine the set of parameters \x t by iteratively setting 

s Rfl (t) = s R (Ut)), 

wherein we identify the parameter x in the iterated integrals with the regularization pa- 
rameter for bare Green functions in QFT and use the representation (|32|) for those Green 
functions. A similar equation is used to determine the set of parameters b t in the iterated 
integrals such that 

r^ b {t) = r R {<f> z {t)). 

If we use a on- shell subtraction at the point z = 1 for example, one simply finds 

A* t = e^(^tog[^)B t (o ; )]). 

Note that this quantity has no essential singularity at x = due to the fact that t'x#^B t 
has the form 1 + 0(x) in our adopted normalization. In this simple on-shell scheme, b t turns 
out to be 

h = mz-*^ x 

and the resulting generalized iterated integral is a well-defined finite quantity which repro- 
duces the renormalized Green function. 

Changing renormalization conditions or schemes is now evidently a convolution of iter- 
ated integrals. For example, transforming to a MS scheme amounts to a convolution similar 
to the one described before, and we would eventually have 

The second Green function on the rhs uses the same fit, h as defined above, the first on 
on the rhs uses the same fi t but determines fif IS using minimal subtracted QFT Green 
functions. Such convolutions were already applied in M and systematically derived in ||. 
They essentially reduce the calculation of renormalized Green functions to a determination 
of appropriate functions B t from which all other results follow via the algebraic structure of 
the Hopf algebra. Note that the convolution also permits us to describe the change of scales 
in bare Green functions; for example 4> z {t) — > 4> pz {t) is described by 

Mt) = [4>z*r Rz (<p pz )]{t), 

which immediately translates to iterated integrals as the convolution T^^i = T^^^by where 
b' t = b t p- #i - t)x . 

In the following we determine some of the B t explicitly, to illustrate the methods. 
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A. Massless Yukawa theory 



Let us apply the previous ideas to massless Yukawa theory, with interaction 

For illustrative purposes we consider a three-loop contribution to the fermion self-energy 
and an analogous one for the vertex part, as sketched in Fig. |5|. 

In order to describe the results, in the context of dimensional regularization, the param- 
eter x is identified with (4 — D)/2 = e and z with the scale factor (J?/p 2 , where p is the 
external momentum and is the renormalization mass scale which arises via the dimen- 
sionality of the bare coupling: g 2 = (— Aix)~ x ix 2x Aira/Y(l — x). The calculations bring in 
supplementary integrals which are readily calculated in any dimension D = 4 — 2x, namely 



g 2 I(p, nx) 



-Wo 



d D k/(2n 



az 



where 



and 



where 



g 2 J(p,nx) = -igl 



Jn+l 



k 2 ((k + p)2(i+™0 4vr (p 2 ) nx ' 

r(l - (n + l)x)r((n + l)x) 
T(l + nx)T(2- (n + 2)x) ' 

r l+p.k/p 2 ]d D k/{2-K) D _ az x j n+1 



k 2 {{k + p)^ 1+nx ) Att {p 2 ) nx ' 

r(2- {n + l)x)T({n + l)x) 
r(l + nx)r(3- {n + 2)x) 



(40a) 



(40b) 



(41a) 



(41b) 



Those functions i n +i, jn+i play the role of the functions B w ^ of the previous considerations. 
Applying the on-shell subtraction which sets z — > 1, Ri, to this graph, yields the following 
result for the self-energy integral: 



a V 
4W 



(z 3x - l)j(j 3 - 2(z 2x - l)j(j 2 + {z x - l)jf 



which is obtained as m[(S , ^ 1 ® id)A(j t (z))} where j t (z) = fl v j w ^ [az plays the role 

of 4> z (t)', A acts on the tree parametrizing j t such that A(j t (z)) = J2jt (1) (z) <E> jt {2) (z) and 
H4jt(z)}-jt0). 

Similarly, the following result for the vertex integral at zero meson momentum is ob- 
tained: 



Gtt) 



(** - l)jfis - 2(z 2x - + (z x - 1)3(1, 



2 / 2 

\i Ip 



in which we evidently see, as before, the (id — R) structure of renormalized Green functions. 
The limit x — > may readily be taken in the renormalized answers, producing 



a 

47T 



IP lP_ u_ 
24 + 16 + 16 



u = Hn 2 / P 2 
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for the self-energy and 

/ ™ \ 3 rrr3 772 rrl 

U = ln(/i 2 /p 2 ) 



u 3 u 2 u 
12 + T + 2" 



for the vertex part. The latter carries an extra weight factor 3 because of the way that the 
interior self-energies can be distributed. 

B. Quantum electrodynamics 

In this case the renormalization scale /i enters through el = (— An)~ x fi 2x Ana/T(l — x), 
where a now has the significance of the fine structure constant. Because we are dealing with 
the massless version, the self-energy assumes the form, = jpA(p 2 ), while the vertex at 
zero-momentum transfer contains two terms (form-factors): 

r>,p) =7^(p 2 )- tn^G(p 2 )/p 2 , 

of which only the first requires (infinite) renormalization. The form factors can be found in 
the usual manner by tracing appropriately: 

Tr[ 7 -T] = 2 2 ~ x [(4-2x)F -xQ] 
Tr[^p.r]=p 2 2 2 -*[.F-S], 

and it is helpful to define the auxiliary integral 

e n (p, nx) - An ^ nx > 

where 

(l-x)r((n + l)a:)r(l-(n + l)a:) 

h n +i = ~ 



2r(2 + nx)r(3 - (n + 2)x) 
We should also exhibit a combination which features prominently in actual computations, 



viz. 



i n+1 + Dh n+1 = T(2 - (n + l)x)r(l + (n + l)x)/r(3 - (n + 2)x)F(2 + nx); D = 4 - 2x. 
Lowest order in a calculations give (z = fi 2 /p 2 ) 

* = (42a) 



,2 / a \ 



,2tt 

leading simply to the on-shell subtracted expressions, 



^=(l-x)^-j^, 

S = / [*i + Dh ^ X i (42b) 
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A 



1]. 



(43a) 



fl,OS 



a 
2^ 



3i 



i»{z x - i; 



2x^V 



(43b) 



It is very easy to continue this process to higher orders for the self-energy function A(p 2 ) 
in the Feynman gauge say; all we need do is to substitute g 2 in massless Yukawa theory 
by 2e 2 in QED, the factor of 2 arising from the gamma- matrix trace. However the vertex 
computations are subtler, because they involve the two form factors T and Q. Here we 
require knowledge of the generic integral, 



Kit 



d u K 



(A7T) D k 2 />+ jk 



fnlti - g, 



ii>+ ft) 

(p + k) 2 



(p + k) 2 ) fc+ Jk 
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a 
2^ 



fn+llv — 9n+l 



(i>+ ft) 
( P + k) 2 



p 2 



(n+l)x 



(44) 



starting with / 
property, 



l,go = 0. In this way we straightforwardly arrive at the recurrence 



v n+ i = M n+ iv n , 

where Vj is a column vector [f il gi) T with v = (1, 0) T and M n is a 2 x 2 matrix with entries 
(M n ) rs given by 



(M n )n 
(M n ) 21 



x) 2 r(nx)T(l + nx) 



r(2 + (n - l)a:)r(3 - (n + l)x) 

r(2 — nx)r(l + nx) 
T(2 + (n - l)x)r(3 - (n + l)x) 



(M n )l2 

(M n ) 22 



r(l — nx)r(na;) 



(45) 



T(l + (n- l)x)r(2 - (n + l)x) 
0. (46) 



We now introduce the matrix-function R with entries Ru = R\ and elsewhere, where 
R\ is again the map which sends z = fi 2 /p 2 — > 1. Further, we use the 2x2 unit matrix 
I such that an insertion of a fermion self-energy amounts to an insertion of matrices of 
the type := I^rjfc into a string of matrices. Then we can express the result for any 
diagram T(t), which iterates one-loop vertices at zero-momentum transfer and also has 
various combinations of one-loop fermion self-energies as subdivergences, by the formula 



r(t) 



n 



D*(«) - 

n w(v) Z 



■w(r)x 



where t is a decorated rooted tree allowing for one of two possible decorations i(v) at each 
vertex: either i(v) evaluates to a one-loop vertex function or to a one-loop fermion self- 
energy. As both decorations are of one-loop order, w(v) is the usual vertex weight discussed 



before. We then set B 



i(v) 
j 



Mj if the vertex v of t is decorated by the one-loop vertex- 
function, and Bj K "' = Aj if the vertex v is decorated by a one-loop fermion self-energy. 
Various examples are drawn in Fig. E. 
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Renormalization now amounts to inserting the matrix R into the string of matrices at 
places which correspond to cuts at the decorated tree. Hence the formulae for the coun- 
terterms and renormalized Green functions remain essentially unchanged. For the Feynman 
graph with the tree t shown in Fig. || we find the bare Green function 

F(t) = AxAiMaVo^" 3 ". 

The counterterm being Zr = — R[r(i) + m{(Zn ® id)(T <g> r)A'(t))], we get 

-RfAxAiMaVo] + 2R[A 1 ]R[A 1 M 2 v ] - RfAijRfAjRfMxVo]. 

Finally the renormalized Green function m[(Zn ® id){T <g> T)A(t)] leads to the on-shell 
subtracted result (aside from a weight factor of 3) 

(^) 3 ft [7, ((M 3 ) n (^ - 1) - 2(M a )ii(* 5ta - 1) + (M 1 ) 11 (z x - 1) 

-(M 3 )ia + 2(M 2 )i2 - (Mi) 12 ) 

- ((M 3 ) 12 z 3x - 2{M 2 ) l2 z 2x + (M x ) 12 z x ) ^ frlv 2 . 

The renormalized contribution may be found by proceeding to the delicate limit x — > 0, 
when the above answer collapses into 



a 
2^ 



7a." 



2U 3 + 3U 2 + 15U-9 
48 



P 



>2U 2 -2U + 3 
16 



u = Mf^/p 2 



Other diagrammatic contributions may be extracted in much the same way; we have 
evaluated many more results for the renormalized self-energy A os and the vertex at zero 
momentum transfer T^ tOS . A useful check on the calculations, that we have carried out, is 
to verify the Ward identity, T{p,p) = dTi(p)/dp fl , since a particular contribution to A will 
produce a combination of vertex graphs with appropriate weights. 



C. A relaxed shuffle product 



Proper iterated integrals satisfy a shuffle identity |I|[TT]. Due to this shuffle identity, one 



can express all tree-iterated integrals representing rooted trees as proper iterated integrals. 
The latter are based on trees without sidebranching and form a closed Hopf subalgebra of 
Ti. The generalized iterated integrals which we used to summarize some restricted results 
of QFT do not obey a shuffle algebra, and hence we get non-trivial representations of the 
full Hopf algebra M. 

However one can easily establish that the leading coefficient of bare Feynman diagrams 
still obeys such an identity, by making use of the expression in terms of generalized iterated 
integrals. For that it suffices to realize that if the \xt become t- independent, we get a proper 
iterated integral. But to leading order, we have 

B > = t^m m > (47) 

with F(t) = l+0(x), by definition; so to leading order, F(t) = 1 Vi, and then B t = l/t'x*^-', 
which can be easily obtained as the iterated integral Gi ;00 (t) defined previously. This explains 
relations known to practitioners of QFT that hold between leading coefficients of the overall 
divergence of graphs representing different trees. 
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VI. COMMENTS AND CONCLUSIONS 



In this and previous papers we have developed a number of mathematical tools which 
clarify and simplify the renormalization procedure for any renormalizable quantum field 
theory. These are encapsulated in formulae such as (11), (13), (15), (26) and (30); we have 
provided two nontrivial examples where these ideas can be fruitfully applied without too 
much effort. Therefore, it is our firm belief and hope that all these tools will enable the prac- 
titioner of QFT to obtain the amplitudes associated with renormalized Feynman diagrams 
by purely combinatorial means. We also anticipate that the connection between these am- 
plitudes and iterated one-dimensional integrals, whose endpoints define the renormalization 
scheme, will also lead to substantial progress in automating the method, since they provide 
a significant step along this direction. The mathematical tools emphasize the fundamen- 
tal connection between geometry, topology, number theory and Feynman diagrams that is 
emerging (T7|jg,|TB]JT^T|. 
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APPENDIX 
A. The Hopf Algebra 

We follow HQ closely. A rooted tree t is a connected and simply-connected set of oriented 
edges and vertices such that there is precisely one distinguished vertex which has no incoming 
edge. This vertex is called the root of t. Further, every edge connects two vertices and the 
fertility f(v) of a vertex v is the number of edges outgoing from v. The trees being simply- 
connected, each vertex apart from the root has a single incoming edge. 

As in [[|, we consider the (commutative) algebra of polynomials over Q in rooted trees; 
hence the multiplication m(t, t') of two rooted trees means drawing them next to each other 
in arbitrary order. Observe that any rooted tree t with root r yields f(r) trees ti, . . ., t/( r ) 
which are the trees attached to r. The unit element of this algebra is 1, corresponding, as 
a rooted tree, to the empty set. 

Let B- be the operator which removes the root r from a tree t: 

B-:t^B-(t)=tit 2 ...tf( r ). (48) 

Fig. |7| depicts an example. Also let B + be the operation which maps a monomial of n rooted 
trees to a new rooted tree t which has a root r with fertility f(r) = n which connects to the 
n roots of ti, . . . , t n . 

B+:t 1 ...t n ^B + (t 1 ...t n )=t. (49) 
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This is clearly the inverse to the action of _B_. For any rooted tree t one has 



B + (B_(t)) = B_(B + (t))=t, (50) 

and Fig. || provides one such example. We further set -B-(ti) = 1, B + (l) = t%. 

We will introduce a Hopf algebra on such rooted trees by taking the opportunity to cut 
such trees in pieces. We start with the most elementary possibility. An elementary cut is a 
cut of a rooted tree at a single chosen edge, as indicated in Fig. |9|. 

Before introducing the coproduct we finally introduce the notion of an admissible cut, 
also called a simple cut. It is any assignment of elementary cuts to a rooted tree t such that 



any path from any vertex of the tree to the root has at most one elementary cut. Fig. |10 
depicts such a situation. An admissible cut C maps a tree to a monomial in trees. If the 
cut C contains n elementary cuts, it induces a map 

<? : f - <7(f) = f[ **■ (51) 

i=l 

Note that precisely one of these trees tj i will contain the root of t. Let us denote this 
distinguished tree by R (t). The monomial which is delivered by the n — 1 other factors is 
denoted by P c (t). The definitions of C, P, R can be extended to monomials of trees in the 
obvious manner, by choosing a cut C l for every tree tj i in the monomial: 



C{t n . 


■ ti ) 

Jn J 


:=C\t n ). 


■ c n (t jn ), 


P c (t n . 


■ti ) 

Jn / 


:=P C \t n ) 


■■■P Cn (t jn 


R (tji ■ 


■■ti ) 

Jn / 


:=R C \t n ) 


■■■R cn (t jn 



Let us now establish the Hopf algebra structure. Following [^,^J we first define the counit 
and the coproduct. The counit e: A — > Q is simple: 

e(X) = 

for any 1/1, 

e(l) = I. 

The coproduct A is defined by the equations 

A(1) = 1®1 (52) 
A(ti ...t n )= A(ti) . . . A(t n ) (53) 
A(t) =t(g)l + (id® B + )[A(B4t))], (54) 

which defines the coproduct on trees with n vertices iteratively through the coproduct on 
trees with a lesser number of vertices. See Fig. 11 for an example. Actually, the coproduct 
can be written as 

A(t) = 1 ®i + t ® 1 + P C (t)(g)R c (t)=:l®t + t®l + A'(t), (55) 

adm. cuts C of t 

which defines A'. 
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Up to this point we have established a bialgebra structure, but it is actually a Hopf 
algebra. Following |2||| we find the antipode S as 

S(l) = 1 (56) 
S(t) = -t- T S[P c {t)}R°(t) = -t-m[(S®id)A'(t)}. (57) 

adm. cuts C of t 

An alternative formula for the antipode, which one may easily derive by induction on the 
number of vertices is 

s(t) = - £ (-ircp c (t)R c ( t ), 

all cuts C of t 

where nc is the number of single cuts in C . This time, we have a non-recursive expression, 
summing over all cuts C and relaxing the restriction to admissible cuts. The overall minus 
sign can be incorporated in the sum if we attach an incoming edge to the root. All cuts 
which remove this edge are then full cuts, all the other ones are normal cuts. A Feynman 
graph corresponding to a tree with n vertices allows then for 2 n cuts. The 2 n_1 full cuts 
deliver the counterterm, the 2 n ~ 1 normal cuts eliminate the subdivergence to deliver the 
result of the i?-bar operation . 

So far we have established a Hopf algebra on rooted trees, using the set of rooted trees, 
the commutative multiplication m for elements of this set, the unit 1 and counit e, the 
coproduct A and antipode S. We call this Hopf algebra 1~Lr. Continuing in the manner of 
we may label the vertices of rooted trees by Feynman graphs without sub divergences, 
in the sense described in the paper and in detail in [|3]{|] . 

Let us also mention that 

m[{3 ® id)A{t)\ = e(t) = = £ S(t {1) )t {2) , (58) 

where we introduced Sweedler's notation A(t) =: J2t(i) ® £(2)> and id is the identity map 
Hr — > TCr. We conclude by defining #(t) to be the number of vertices of a rooted tree t. 
This extends to a monomial of rooted trees in the obvious manner: = J2i 



B. The Hopf Algebra Structure of Graphs and Forest Formula 

The results of [|J show that for each Feynman graph T, we obtain a sum of associated 
rooted tree T r and a coproduct given by 

A(T r ) = 1 ® T r + T r <g> 1 + £ T 7 ® T r/7 . (59) 

Here, Tp is a sum of rooted trees with decorations, primitive elements in the Hopf algebra 
of rooted trees, which are obtained from Feynman graphs without sub divergences. As the 
map Tp <-> T is one-to-one, we can directly formulate the Hopf algebra on graphs T, as in 
section II. 

To the coproduct ( |59| ) belongs an antipode given by 

S(T r ) = -T r -J2 5[T 7 ]T r/7 , (60) 
7cr 
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as one straightforwardly checks. Because it is an antipode in a Hopf algebra of rooted trees, 
it can be written as a sum over all cuts. Set T r = J2i Ti for some decorated rooted trees T. 
Then, 



S(T r ) = £ E (-l) nc ^' 



(Ti)R Ci (Ti). 



(61) 



i all cuts d of Ti 



Each such cut corresponds to a renormalization forest, obtained by boxing the corresponding 
subgraphs in T, and vice versa 0. 

Now, let be a Q-linear map which assigns to Tr the corresponding Feynman integral. 
Further, let 4>r ; = tr o be a map which assigns to Tr the corresponding Feynman integral, 
evaluated under some renormalization condition R. Hence, from Tr we obtain via <fi a 
Feynman integral 0(Tp) in need of renormalization. tr modifies this Feynman integral, 
in such a way that the result contains the divergent part of this integral. Essentially, tr 
extracts the divergences of 4>(Tr) in a meaningful way |I2[ . Hence, as tr isolates divergences 
faithfully, differences (id — r/j)(0(Tr)) eliminate infinities in Feynman integrals. Depending 
on the chosen renormalization scheme R, one can adjust finite parts to fulfill renormalization 
conditions. A detailed study of this freedom from the Hopf algebra viewpoint can be found 



We remind the reader of Sweedler's notation: A(T r ) = I^Tpm <S> ^r(2)- Let us consider 
the antipode e(T r ) in the same notation: 



in 



= e(T r ) = ££(T r(1) )T r(2) . 



The above map vanishes identically, and it can also be written as 



m[(S <g> id)A(Tr)] = e(T r ) = 0. 



But this map gives rise to a much more interesting map, by composition with 0, 



T r - := m[{S R ® zd)(0 ® 0)A(T r )], 



for it associates the renormalized Feynman integral Tr 0,0] to the Feynman graph T repre- 
sented by a unique sum of rooted trees. 
Its usual definition, 




(62) 



7cr 



is recovered if we define 




(63) 



7'C7 



This map is derived from the antipode 



(64) 



7'C7 
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Using to lift this to Feynman graphs, and using the freedom to alter corresponding analytic 
expressions according to renormalization schemes R one ends up with (|63l). 
Note that if one defines 



57 



s. 



one has Sr o 



S and hence 



5fl[0(T 7 )] = t r 



-0(T 7 ) - 2 ^(5[T y ])0(T 7/y 

y C7 



(65) 



Hence, in accordance with [@-|]] we determine the Z-factor of a graph 7 as derived from 
the antipode in the Hopf algebra of rooted trees. In (|62|) , we recovered the original forest 
formula in its recursive form. The non-recursive form is recovered with the same ease, using 
(0) instead of 
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FIGURES 



FIG. 1. A Feynman graph, its divergent subgraphs, its forests, and the corresponding tree with 
appropriate full and normal cuts. ^From this one calculates the antipode and the renormalized 
Green function. 

FIG. 2. Forests and counterterms as cuts and antipodes. Each forests (thick grey lines) contains 
an analytic expression which has to be evaluated using R. 

FIG. 3. Vertex weights and tree factorials. For the tree t given at the left in the upper row 
we indicate the vertex weights attached to each vertex. We also give the Feynman diagram corre- 
sponding to a chosen common one-loop self-energy decoration at each vertex. 

FIG. 4. Generalized vertex weights and tree factorials which can appear if there are decorations 
of arbitrary loop order. The factorial of this decorated tree is 2 x 5 = 10. 

FIG. 5. A Yukawa self-energy at three loops. 

FIG. 6. Examples of Feynman graphs, their decorated rooted trees and the corresponding 
matrix calculus. From top to bottom, the bare function T(t), evaluates for the indicated decorated 
rooted trees to Miv z~ x , AiM 2 v z~ 2x , AiM 2 M 3 vo-z~ 3:E . 

FIG. 7. The action of B on a rooted tree. 

FIG. 8. The action of B + on a monomial of trees. 

FIG. 9. An elementary cut c splits a rooted tree t into two components, the fall-down P c {t) 
and the piece which is still connected to the root, R c (t). 

FIG. 10. An admissible cut C acting on a tree t, leading to a monomial of trees. One of the 
factors, R c (t), contains the root of t. 

FIG. 11. The coproduct, worked out for the trees t\, t 2 , ts 1: ts 2 , from top to bottom. In the last 
line we give one full and all normal cuts of the tree. 
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FIG. 3. 
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FIG. 6. 
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